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ABSTRACT 


We perform a series of 2D smoothed particle hydrodynamics (SPH) simulations 
of gas accretion onto binaries via a circumbinary disc, for a range of gas temperatures 
and binary mass ratios (q). We show that increasing the gas temperature increases the 
accretion rate onto the primary for all values of the binary mass ratio: for example, 
for q = 0.1 and a fixed binary separation, an increase of normalised sound speed 
by a factor of 5 (from our “cold” to “hot” simulations) changes the fraction of the 
accreted gas that flows on to the primary from 10% to ~ 40%. We present a simple 
parametrisation for the average accretion rate of each binary component accurate 
to within a few percent and argue that this parameterisation (rather than those in 
the literature based on warmer simulations) is relevant to supermassive black hole 
accretion and all but the widest stellar binaries. We present trajectories for the growth 
of q during circumbinary disc accretion and argue that the period distribution of stellar 
“twin” binaries is strong evidence for the importance of circumbinary accretion. We 
also show that our parametrisation of binary accretion increases the minimum mass 
ratio needed for spin alignment of supermassive black holes to q ~ 0.4, with potentially 
important implications for the magnitude of velocity kicks acquired during black hole 
mergers. 


1 INTRODUCTION 

Understanding the accretion of material onto binaries is 
important for many astrophysical problems. Super massive 
black hole binaries are thought to occur frequently as a result 
of large galactic mergers. However, observationally detect¬ 
ing black hole binaries remai ns a challenge (|Artymowicz fc 


Lubow|1996||Dotti et al.|2012[|Schnittman|20l3 i Bogda novic 

201 5| . As such, understanding the electromagnetic signa¬ 


tures of gas flows onto binaries is important for improving 
detection . The relative accretion rate of gas by the two bi¬ 
nary components is also important in setting the alignment 


of black hole spins ( 

Bogdanovic et al. 2007} Lodato & Gerosa 

2013 Gerosa et al. 

2015). The alignment of black hole spins 


sets the recoil velocity of the post-merger black hole and can 


influence the host galaxy through feedback processes (Blecha 


& Loeb 2008). Each black hole’s spin is aligned with its 


disc by the Bardeen-Petterson effect (Bardeen & Petterson 


1975 ), on a time scale which depends upon the accretion rate 
through the disc ( Gerosa et al.|2015| . As such, understand¬ 
ing the relative strength of primary (Mi) and secondary 
(M 2 ) accretion is vital for determining the likelihood of spin 
alignment. 

The distribution of stellar mass within binaries is also 
sensitive to the details of binary accretion. Given that a 
large fraction of stars are in binaries, this problem is of great 
significance ( Close et al.|2003| Basri &; Reiners|2006| Fischer 


&; Marcy||1992| |Duquennoy &; Mayor||1991| |Preibisch et al.| 


1999; |Mason et ah|1998| . It is well known that the initial 
mass of a protobinary is only a fraction of the mass of the 


parent cloud (except for the widest binaries) (Boss 1986 
Bonnell & B ate|199 4). As most of the mass not in the initial 


protobinary is of greater specific angular momentum, the 
protobinary can only continue to grow by accreting material 
through a circumbinary disc. The ultimate mass ratio of the 
binary (q = M 2 /Mi, where subscripts 1 and 2 refer to the 
primary and secondary, respectively) is set by the accretion 
rates onto the primary (Mi) and secondary (M 2 ) through 
the circumbinary disc. 

It is possible that the circumbinary disc that feeds ma¬ 
terial onto the binaries may be misaligned with the binary 
orbit. Accretion via a misaligned circumbinary disc does not 
change the component of the binary that accretes the most 
material (Nixon et al. IpOTfl [UTS] |Dunhill et aTpoH) , al- 
though the accretion rate is potentially higher for a mis¬ 
aligned circumbinary disc ( Nixon et al.|2013 ). 

Cluster level simulations have been used to predict 
properties of binaries (Goodwin et al. 2004| Offner et al. 


2009 Bate|200~9 Delgado-Donate &; Clarke|2008 ). They find 

mass distributions peaked towards high values of q , implying 
preferential accretion onto the secondary. Dedicated studies 
of isolated binaries, which allow for much greater resolution, 
have also been performed using both smoothed particle hy- 
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2 Young and Clarke 


drodynamics (SPH) in 3D ( Bate &; Bonnell|1997| Dotti et al. 
2010 ) and 2D (|Young et al.||2015 ), and using grid codes i 


2D 


(2015) and inject material at a fixed boundary far from the 


Han awa et al.||2010 de Val-Borro et al. 2011 Farris 


et al. 2014). The overwhelming majority of these studies 


find a strong tendency for the majority of the mass to flow 
onto the secondary. This is because material flowing from 
the circumbinary disc encounters the secondary first, since 
its orbit is further from the centre of mass of the system. In¬ 
deed, the same preference for accretion onto the secondary 
can be found in simulations of ballistic (non-collisional) ma¬ 
terial with specific angular momentum greater than the bi¬ 
nary (B ate||199 7). 

However, simulations have also shown that gas temper¬ 
ature can moderate this preference for accretion onto the 


secondary (Dotti et al. 2010 Young et al. 2015). This is 


because at higher temperatures, material travels from the 
circumbinary disc onto the binary along a wider range of 
trajectories, allowing some of it to orbit the edge of the sec¬ 
ondary’s Roche lobe and accrete onto the primary (Young 


et al. 2015). Understanding the dependence of relative ac¬ 


cretion rate on both gas temperature and binary mass ratio 
is vital for providing accurate predictions of black hole spin 
alignment and the distribution of q in stellar clusters. Un¬ 
derstanding the temperature dependence of accretion rates 
is particularly important for black hole applications, where 
the ratio of sound speed to binary orbital speed is too low 
to be modeled by current simulations ( Gerosa et al.||2015 ). 

In this paper we present a series of SPH calculations of 
the accretion of gas onto binaries at a range of gas temper¬ 
atures and mass ratio. We use our simulations to propose a 
simple parametrisation of the relative accretion rate of bina¬ 
ries, as a function of both gas temperature and mass ratio. 
In Section [2] we describe the details of our model and the 
computational set up. In Section [ 3 ] we present the results of 
our simulations. In Section [4] we discuss the interpretation 
of our results, their implications. Finally we provide some 
concluding remarks in Section [5] 


2 MODEL DETAILS 

We model a binary with mass ratio q = M 2 /Mi, separation 
a, total mass M = Mi + M 2 on a anti-clockwise circular 
orbit. All simulations are performed in 2D using a modified 
version of the SPH code GADGET2 (Springel 2005). We 


choose to use 2D rather than 3D for comparison to earlier 
work (Young et al. 2015) and to reduce computational costs 
(for N particles, resolution scales as iV -1//2 in 2D and iV -1//3 
in 3D). 

Physically, accretion onto a binary is mediated by a cir¬ 
cumbinary accretion disc. As viscosity in this disc causes it 
to spread inwards, material flows from the edge of the cir¬ 
cumbinary disc onto the binary. Self consistent modelling of 
this viscous spreading, which is primarily driven by an ef¬ 
fective viscosity resulting from turbulent processes (Pringle 


1981 Shakura & Sunyaev 1973), is computationally chal¬ 


lenging for binary accretion modeling and has not to date 
been extended to studying accretion onto the individual bi¬ 
nary components (Shi et al.|2012 ). Furthermore, the details 
of the accretion in the circumbinary disc are unlikely to af¬ 
fect the steady-state accretion rate of the binary. Given this, 


binary, R — Ri n fd , with fixed specific angular momentum, 
j = jinfCGMa. The initial radial velocity of material at the 
outer edge is set such that the material is initially marginally 
gravitationally bound. This allows us to provide a continu¬ 
ous supply of material to the binary, without the computa¬ 
tional expense of simulating a large circumbinary disc. 

The mass of the gas is taken to be negligible compared 
to the binary and so the self-gravity of the gas is ignored. For 
the same reason, the binary masses and orbits remain fixed 
throughout our simulations. Any particle which is within 
R — Raced of a binary particle is removed from the simula¬ 
tion and the number of removed particles per binary particle 
is recorded (we set R acc — -02 in all our simulation). The 
gas temperature is modeled using an isothermal equation of 
state, which is parametrised in terms of c: 


GM 

C = C W — 


(1) 


c is then the sound speed normalised to the orbital speed 
of the binary (it is also H/R of the circumbinary disc at a 
radius of a). 

For all simulations performed here we set ji n f = 1.2 
and Rinf = 20 which ensures a constant supply of material 
onto the binary with specific angular momentum similar to 
what would be expected for material at the inner edge of a 
circumbinary disc (Ba te &; Bonnell|T9 97 j Bate 2000). Lower 
values of ji n f would be unrealistic for accretion via a disc, 
while higher values of ji n f would cause a disc to form at 
a radii where material cannot flow onto the binary until 
viscous spreading moved the boundary inwards. We inject 
particles at the boundary at a rate of N = 500 per dynamical 
time, which ensures the edge of the circumstellar discs are 


well resolved (Young et al. 2015). As the temperature of 


the gas has been shown to affect the relative accretion rate 
of binary components, we explore values in the range c = 
0.05 — 0.25. As lower values of c produce thinner discs which 
require more resolution elements to resolve, we were unable 
to probe values of c < 0.05. We run simulations with 0.1 < 
q < 1.0 to explore the effect of q on relative accretion rates. 
Smaller values of q < 0.1 were deliberately excluded as they 
are more computationally demanding and the nature of the 
accretion flow has been shown to change qualitatively in this 
range. That is, for q < 0.1, the accretion transitions from 
gas flowing onto a binary in a low gas density cavity, to gas 
flowing through a large circumprimary disc, with a large 
embedded body clearing a gap (D’Orazio et al. 2013). Each 
simulation was run for at least 500 binary dynamical times, 


we choose to follow Bate & Bonnell (1997) &; Young et al. 


by which time they have settled into a steady-state (Young 
et al.|2015 l. 

All analysis was performed in a frame corotating with 
the binary, translated so that the primary is at the origin 
and the secondary is at (x,y)=(-l,0). In this reference frame 
it is possible to define a modified potential, commonly called 
the “Roche potential”, as: 

_ GMi _ GMz 1 2 

|r — ri| |r — r 2 | 2 V ; V ; 

where ft is the angular momentum vector of the binary, r 
is the position vector and subscripts 1 and 2 refer to the 
primary and secondary respectively. Surfaces of equipoten- 
tial of 4> and Lagrange points (points of zero gradient) are 
shown in Figure [Tj 
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Binary accretion rates: dependence on temperature and mass-ratio 3 


Overview of Roche lobes and Lagrange points 



Figure 1 . Equipotential surfaces of <£>, the modified potential 
(gravitational potential plus term to include centrifugal force of 
non-inertial frame) for the values of <f> obtained at the three La¬ 
grange points L1,L2 and L3. The colour map shows <f> elsewhere. 
The stars are marked in red (secondary) and green (primary) and 
the centre of mass of the system is marked in blue. The green and 
blue lines are the locations where the x and y components of the 
force resulting from this potential (i.e., — V<f>) are zero. 


ID 

q 

c 

^end 

1 

0.1 

0.05 

720 

2 

0.2 

0.05 

500 

3 

0.3 

0.05 

640 

4 

0.4 

0.05 

630 

5 

0.5 

0.05 

690 

6 

0.6 

0.05 

720 

7 

0.7 

0.05 

730 

8 

0.8 

0.05 

840 

9 

0.9 

0.05 

800 

10 

0.1 

0.25 

2000 

11 

0.2 

0.25 

2000 

12 

0.3 

0.25 

2000 

13 

0.4 

0.25 

2000 

14 

0.5 

0.25 

2000 

15 

0.6 

0.25 

2000 

16 

0.7 

0.25 

2000 

17 

0.8 

0.25 

2000 

18 

0.9 

0.25 

2000 

19 

0.1 

0.10 

542 

20 

0.1 

0.15 

521 

21 

0.1 

0.20 

515 


Table 1 . Parameters used for simulations in this paper. For each 
simulation particles were injected at a rate of roughly N « 500 per 
dynamical time (y/a 3 /GM), the accretion radius was set to 0.02a, 
particles were injected from 20a with specific angular momentum 
1.2VGMa. 


Table ^ lists the parameters used for all simulations in 
this paper. Details of the exact version of the SPH code 
used along with initial conditions and parameter files can 
be found in Section [6] 



Figure 2. Steady-state values of A = M\/M for simulations with 
cold gas (c = 0.05), for a range of values of q. The central line 
shows a moving average over 3 binary orbits of A while the shaded 
region shows the one standard deviation range on A within this 
range. The first 200 binary dynamical times are omitted to ensure 
that only steady-state accretion is shown. 


3 RESULTS 

All our simulations undergo a transient settling phase, which 
is dependent on the details of the initial conditions, before 
settling into a steady-state after ~ 30 binary orbits. To 
exclude the physically uninteresting effects of this settling 
phase, we excluded the first 200 dynamical times (where 2n 
dynamical times is one binary orbit) from all the results that 
follow. 

The evolution of the binary mass ratio q and the rel¬ 
ative strength of primary/secondary accretion flows can be 
characterised using the parameter, 

r -i / s- e ^ <A -•*»■> (3) 

which is the fractional change in q per fractional change in 
M. This can be re-written as, 

r= ^ 7^ (1_a(1 + 9)) (4) 

where A = Mi/M and M — Mi + M 2 . In our simulations 
we measure Mi and M 2 by counting the number of particles 
accreted onto the primary and the secondary. 

The evolution of A as a function of time and q is shown 
for simulations with cold gas (c = 0.05) in Figure [ 2 ] A is 
smoothed using a moving average over three binary orbits, 
which is roughly the orbital period of the edge of the binary 
cavity. 

Figure [3] shows the same information for the hot sim¬ 
ulations (where c = 0.25). Unlike the cold simulations the 
accretion rate A is significantly more variable, making it dif¬ 
ficult to discern any trend in the average accretion rate with 

q- 

Figure [4] show the distribution of A in the steady-state 
for all simulations. This allows the trends in the average 
accretion rate with q and c to be more clearly visualised. The 
cross in each distribution shows an estimate of the average 
steady-state A, given by f Mi/ f M. 

Figure [ 3 ] also shows the distribution of steady-state ac¬ 
cretion rates for a series of simulations with q — 0.1 and 
gas temperatures intermediate between hot (c = .25) and 
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Figure 3. The same as Figure [2] but for simulations with hot gas 
(c — 0.25). 



Figure 4. The distribution of steady state A for all simulations. 
Each simulation is represented by a violin plot, where more com¬ 
mon values of A correspond to a greater horizontal extent. The 
crosses show the average steady-state accretion rate across the 
entire simulation (excluding t < 200 td yn )i f Mi/ f M. Hot simu¬ 
lations are shown in red, cold simulations in blue and intermedi¬ 
ate simulations in yellow. The blue line shows the parametrisation 
A = q/2 and the green line shows A = —|-j- 


cold (c = .05). These additional simulations show how the 
accretion rate changes with gas temperature in more detail. 


4 DISCUSSION 


4.1 Variability of A 


All simulations show a significant amount of temporal vari¬ 
ability in A. This variability is driven by fluctuations in the 
amount of material flowing into the binary cavity, which has 
been identified in many previous simulation of gas flow onto 
binaries ( D’Orazio et al.|2013||Qchi et al.|2005| Young et al. 


2015 [Artymowicz &; Lubow 1996; MacFadye n &; Milosavl- 


jevic 2008 Shi et al.||2012). 

The physical mechanism for the variability in supply of 
material onto the binary was explained in detail b y |Mac-| 
Fadyen &, Milosavljevic (2008) &; |D’Orazio et alT (2013). 


They found that the edge of the circumbinary disc was ec¬ 
centric (i.e., the cavity wall is not circular). Because of this, 
the amount of material flowing inwards fluctuates on the 
time scale of, 1) the binary orbit and 2) the precession time 
scale of the disc edge. 

Because changes in the rate of flow onto the binary are 


Figure 5. The vertical axis is the standard deviation in the rate 
of inflow onto either Roche lobe, divided by the particle injection 
rate at the outer boundary (i.e., the standard deviation of the 
inflow rate into the cavity). As with all other quantities, the first 
200 dynamical times were excluded from the calculation so that 
only the steady state is probed. 


driven by disc eccentricity, the magnitude of the variability 
(and the importance of the different time scales) varies with 
q. Generally, a larger q is able to excite a more eccentric 
cavity edge, which leads to greater variability (see |D’Orazio 


et al. (2013) for details of the physical mechanism). Fig¬ 
ure [5] shows the variability in the amount of gas flowing in¬ 
wards past the Roche equipotential surface passing through 
the LI point. Although there is significant variability in the 
infall rate, the amount of variability does not change sys¬ 
tematically with q. This is likely a limitation of our deci¬ 
sion to inject material with constant angular momentum 
j = 1.2V GMa , which prevents the disc edge reaching high 


eccentricities (MacFadyen & Milosavljevic (2008); Shi et al. 
(2012) found e ~ 0.1 —0.4 for q— 1.0) by providing a flow of 
material that will push the disc towards a circular flow be¬ 
yond the circularisation radius of R — j% n fd. On the other 
hand, a trend towards increased variability at higher gas 
temperatures can clearly be seen. 

Unlike Figure [5] Figure [4] shows a clear trend towards 
the accretion rate onto either component of the binary be¬ 
coming less variable as q increases. This is because A (and 
the accretion onto the secondary, 1 — A) measures the ac¬ 
cretion rate through the inner edge of the discs around each 
component of the binary. Material flows from the edge of the 
circumbinary disc onto either the circumprimary or circum- 
secondary disc. To reach the inner edge and be accreted, it 
must be “processed” by the circumstellar disc. That is, ma¬ 
terial that joins the circumprimary disc is only registered 
at the inner edge once that material has been moved there 
by the circumprimary disc’s viscosity. For discs with lower 
viscosity (e.g. colder discs) or with a greater spatial extent 
(e.g. discs with higher q), variability in the supply of mate¬ 
rial at the outer edge will be more smeared out by the time 
the inner edge is reached. 

Given this, the variability in Figure is likely driven 
largely by numerics, with simulations with hotter gas and/or 
larger discs (i.e. higher q) transmitting more of the vari¬ 
ability in supply of material striking the discs around each 
binary component through to the inner edge where it is 
accreted. Physically realistic discs, with significantly lower 
viscosity than the numerical viscosity of our simulations, 
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greater dynamic range in disc radius and lower gas tem¬ 
peratures, will have much less variability in the accretion 
rate onto the components of the binary than our simula¬ 
tions. However, the average accretion rate over many binary 
orbits should be independent of the amount of numerical 
variability. This can be seen in the simulations of [Young| 


et al. (2015), where the magnitude of the variability de¬ 


creases with resolution but the average accretion rate is un¬ 
changed. 


4.2 Parametrisation of A 


The most important quantity to understand for our purposes 
is the steady-state value of A, which uniquely determines the 
fractional accretion rate T via Equation^] For the reasons 
discussed in the previous section, it is the steady state av¬ 
erage of A which is important in determining the physical 
properties of the binary. 


It was recently suggested by Gerosa et al. (2015), that 
A can be parametrised as, 


A = 


4 + 1 


( 5 ) 


based upon the data of Farris et al. (2014), who simulated 


accretion onto binaries using a grid-based code for a range 
of values of q and intermediate temperature gas (c = 0.10). 

The green line in Figure [4] shows this parametrisation. 
It is clear that this expression is appropriate for binaries 
accreting gas of intermediate temperature (understandable, 
given the source data), but provides a poor estimate for the 
accretion of hot or cold gas. 

By symmetry, any parametrisation must have A = 0.5 
when q — 1. Based on the data in Figure [4] we propose the 
simple parametrisation 


A = - + m(q - 1) 


( 6 ) 


where m is a function of c and m ^ 0.5 as c ^ 0. The cold 
limit, where A = q/ 2, is shown in blue in Figure [4] and is 
accurate for the cold simulations to within a few percent. 

To asses how m in Equation [6] varies with c, we calcu¬ 
lated the best fit value of m to the data using least squares 
minimisation, which we plot in Figure [6] It is clear that an 
asymptotic limit of m = 0.5 is appropriate for the cold limit 
and that m decreases as a strong power of c. We have avoided 
giving a functional form for m given that we have relatively 
few data points in c and the values of m for c = 0.10,0.15 
and 0.20 are fit using only one sampling in q. 

Moreover, it is the cold limit which is most applicable 
to the bulk of astrophysical systems. In the case of protobi¬ 
naries, c is likely to lie in the range ~ 0.01 — 0.1 (Young et al. 


2015). For black hole binaries, normalised sound speeds as 
low as .001 are possible ( Gerosa et al.|2015 ). In either case, 
our results suggest that the binary accretion rate should be 
well approximated by, 


r - 


(i — < ? 2 )(2 + g) 

2 q 


( 7 ) 


As well as being in good agreement with the results of Sec¬ 
tion [3] this parametrisation is also consistent with the find¬ 
ing of earlier simulations of cold gas accreting onto binaries 
( Bate &; Bonnell||1997 Young et al.|2015 ). 
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Figure 6. The best fit value of the parameter m in the parametric 
model for A given by Equation [6] For each value of c Equation [6] 
was fit using the average steady state value of A for each q. Note 
that the values for c = 0.10,0.15 & 0.20 are fit using just one data 
point. 


4.3 Implications 


For all values of parameter space that we explored, we find 
T > 0, implying that accreting material onto a binary al¬ 
ways brings its mass ratio closer to 1. It has been shown 
that when q < 0.1, a region of parameter space we chose 
not to explore, the physical behaviour of the accretion flow 
transitions to that of a planet embedded in a disc, caus¬ 
ing T to fall ( D’Orazio et al.|[2013 Farris et al.||2014 ). It is 
possible that for very small values of g, T may become neg¬ 
ative, but this would need to be confirmed with specialised 
simulations. 

By integrating Equation [7] we can determine how q 
evolves as a function of accreted mass. This is shown for 
both hot and cold simulations (dashed and solid lines, re¬ 
spectively) in Figure [7] Note that this plot is somewhat dif¬ 
ferent from those presented in |Bate| (2000), in which a proto¬ 
binary was evolved subject to infall from a variety of parent 
cores. The resulting evolution in q was then sensitive not 
only to the density and angular velocity profile of the cores, 
but also to the assumptions made about protobinary migra¬ 
tion. Specifically, Mi/M 2 depends on the specific angular 
momentum of the infalling gas relative to the binary , which 
changes as the binary migrates. Here we define the initial 
binary as corresponding to the point at which it begins to 
accrete from a circumbinary disc. From this point on, the 
material accreted onto the binary must always be of similar 
specific angular momentum to the binary and so q simply 
evolves along the appropriate path in Figure [7| Figure [7] is 
therefore readily applicable to either the stellar binary or 
black hole case. 


4.3.1 Application to stellar binaries 


Main sequence binary star surveys of solar type stars reveal 
a remarkable populati on of “twin-like” binaries (q > 0.95; 


)puiation or twm-nke binaries {q > 
(2010]))[^] In the light of Figure[7]and 


Raghavan et al.| (2U10p)| x | In the light of Figure[7]and given 
the lack of obvious mechanisms that otherwise produce such 


1 The existence of these “twin-like” binaries is not in doubt, how¬ 
ever, there may be some under-representation of the lowest q bi¬ 
naries as their abundance must be corrected in surveys which are 
incomplete (Go ldberg et al.| 20'03|). 
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Figure 7. Change in q as a function of accreted mass over initial 
binary mass, for cold gas with A = q/2 (solid line) and hot gas 
with A = 0.17g + 0.33 (dashed line). 


strongly correlated masses, it is natural to ascribe these ob¬ 
jects to mass equalisation by circumbinary disc accretion. 
Figure [7] implies that such twin-like systems can be pro¬ 
duced from a wide range of initial q values, provided that 
the binary accretes of order several times its own mass from 
a circumbinary disc. On the other hand studies of binary mi¬ 
gration during circumbinary disc accretion all concur that 
the (inward) migration time is less than or of the same or¬ 
der as the mass doubling time (MacFadyen & Milosavljevic 


2008) |Shi et al .|2012 |Cua dra et al. 2009) . This immediately 

implies that “twin-like” systems have migrated inwards by 
a large factor. A further corollary is that the parent sys¬ 
tem (binary plus circumbinary disc) would in the past have 
contained an angular momentum that is considerably larger 
than that contained in the present day binar }0 

The period distribution of “twin-like” binaries amply 
supports this picture. The angular momentum budget of 
pre-stellar cores ( Goodman et al.||1993 ) is only sufficient to 
form binaries out to separations of order 10 3 A.U. [[] Given 
that the angular momentum in the final binary is only a 
fraction of the angular momentum needed to drive a binary 
into a twin-like state, it is unsurprising that the maximum 
separation of twin-like objects is considerably less than this 
(~ 100 A.U.). On the other hand, low q systems are rela¬ 
tively uncommon at small separations (< 5 A.U.); this again 
is consistent with a picture in which the bulk of close bina¬ 
ries are driven to small separation by circumbinary accretion 
(rather than simply arising from cores with exceptionally low 
specific angular momentum). 

Figure [7] also shows that more accreted mass is needed 
to equalise the mass of binaries when the gas is “hot”. As 
the normalised sound speed and seed protobinary mass as 


2 The simulations presented here have been artificially set up 
so that the specific angular momentum of the infalling material 
only modestly exceeds that of the binary in order to mimic the 
inner parts of a circumbinary disc in a computationally efficient 
manner; in reality, however, the inner regions of such a disc are 
fed by viscous redistribution within a reservoir of material whose 
mean specific angular momentum far exceeds that of the binary 
and this becomes increasingly true as the binary migrates inwards. 

3 Wider binaries are attributed either to cluster mediated capture 
( |Kouwenhoven et al.|2010| |Moeckel &; Clarke|2011| ) or to orbital 
reconfiguration in multiples producing wide companions on high 
eccentricity, low angular momentum orbits (|Reipurth & Mikkola| 
2012). 


a fraction of total core mass both increase with binary sep¬ 
aration ( Young et al.||2015| Boss||1986| |Bate||2000| , binaries 
evolve less far to the right (as less Mf/Mi gas mass is avail¬ 
able to be accreted) and along slower tracks (as the gas 
is “hotter”) in Figure [7] As such, we expect small separa¬ 
tion binaries to have q distributions strongly biased towards 
< 7 = 1 , while wide binaries should show a much weaker pref¬ 
erence for equal mass binaries, again consistent with the 
picture described above. 


4-3.2 Black hole binaries 


Gerosa et al. (2015) examined whether circumbinary accre¬ 


tion onto black holes can align their spins with the circumbi¬ 
nary disc on the timescale on which they are driven inwards 
by circumbinary disc torques. This analysis involved pre¬ 
scriptions for total binary accretion rate, disc mediated mi¬ 
gration, accretion induced spin alignment and the relative 
accretion rates onto the binary components as a function of 
q. The latter prescription was derived from the simulations 


of Farris et al. (2014). In conjunction with the other model 


assumptions it was concluded that whereas the secondary 
always aligns with the circumbinary disc, the primary re¬ 
mains misaligned for q < 0.2. Misaligned spins at the point 
of binary merger have important consequences for the recoil 
velocities resulting from the merger event. 

If we repeat this analysis but instead use our prescrip¬ 
tions for differential mass acquisition (which are more ap¬ 
propriate for the ‘cold’ conditions expected in the black hole 
case) we find that it is easier for black hole spins to remain 
misaligned: for the fiducial model this increases the maxi¬ 
mum q value for misalignment from q = 0.2 to q = 0.4. This 
is because we find a stronger preference for accretion on to 


the secondary than in the warmer Farris et al. (2014) sim¬ 


ulations; consequently the primary is relatively starved of 
accretion and can remain misaligned at higher q values. Our 
result thus strengthens the result of Gerosa et al. (2015), 


implying a significant population of black holes whose spins 
are expected to be misaligned at the point of merger. 

When black holes with misaligned spins merge, the 
merger product receives a velocity “ kick”, the maximu m 
value of which scales as 77 = < 7/(1 + q) 2 (Lousto & Zlochower 
2013). For binaries with q > 0.4 the maximum merger kick 


velocity may be comparable with the escape velocity of typ¬ 


ical elliptical galaxies (Merritt et al. 2004). Moreover, the 


fact that binaries may enter the gravitational wave inspi¬ 
ral regime with misaligned spins and relatively large q has 
potential implications for the richness of the gravitational 
wave signature potentially detectable by e-LISA. 

There is however a caveat to this. The prescriptions of 


Gerosa et al. (2015) imply that the binary mass doubling 


timescale is much greater than its migration timescale and 
in this limit it is appropriate to consider spin alignment at 
constant q. However, recent work employing ‘live’ MRI tur¬ 
bulence in the gas flow finds a mass doubling time which 


is only modestly in excess of the migration time (Shi et al. 


2012). If this is correct then migration by a large factor (i.e. 
the ~ 10 fold shrinkage that is required before migration 
due to gravitational wave emission becomes effective), also 
implies a significant change in binary mass (and thus, by 
Figure [7] the value of q). Indeed we have seen above that 
in the stellar case there is ample observational evidence for 
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a preference for high q pairs at small separation. It is not 
clear in the black hole case whether there is a sufficiently 
large mass reservoir in the circumbinary disc to effect a sig¬ 
nificant change in q and it is also unclear whether (for a 
restricted mass reservoir) the relatively strong accretion on 
to the binary limits how far it will migrate. Pending further 
theoretical work in this area, we cannot rule out the pos¬ 
sibility that even low q pairs might not eventually attain a 
high enough value of q for the primary to also align with the 
circumbinary disc. 


5 CONCLUSIONS 


In this paper we have shown that the mass accretion rates of 
the primary and secondary components of binaries depend 
on the gas temperature. We have shown that for all values 
of q , increasing the gas temperature allows more material 
to accrete onto the primary component of the binary. The 
relative accretion by the primary, Mi/M, is also found to 
increase at least quadratically with the sound speed of the 
accreted gas. 

Although our simulations show a range of levels of vari¬ 
ability of accretion onto the star, we argue that it is only 
the higher than physical disc viscosity and lower than phys¬ 
ical dynamic range in disc radii that allows variability in 
the mass influx at the Roche lobe to be manifest as vari¬ 
able accretion onto the star. For realistic physical discs any 
variability in the inflow rate will be smoothed out by the 
circumprimary/secondary disc and the variability in the ac¬ 
cretion rate onto the binary components will be low. As 
such, it is the average steady-state accretion rate that is the 
physically relevant quantity in our simulations. 

We have proposed a simple parametrisation of the rel¬ 
ative accretion rates onto the two components of the binary 
and shown that it reduces to, 


Mi _ q 
M ~ 2 

and 

M 2 _ 2 — q 
~M 2 


( 8 ) 

(9) 


when the accreted gas is cold. 

We argue that the period distribution of “twin-like” 
stellar binaries provides strong evidence that these objects 
acquired the bulk of their mass through circumbinary accre¬ 
tion and that they have migrated inwards by a large factor 
from their initial birthplaces. Applying our parametrisation 
to super massive black hole binaries with existing models 
for spin alignment we find that black hole spins’ remain un¬ 
aligned when q is less than 0.4. However, this value may be 
decreased if the mass accreted during the black hole merger 
is sufficient to drive significant evolution in q. 


6 MATERIALS & METHODS 

In the interests of reproducibility and transparency, 
all the code needed to reproduce this work has been 
made freely available online at https: //bitbucket. org/ 
constantAmateur/binaryaccretion See the readme file in 
this repository for further details. 


All figures in this paper were generated using the 
python package matplotilb (Hunter 2007). 
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